Improved description of 34 ' 36 ' 46 Ar(p,d) transfer reactions 
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An improved description of single neutron stripping from 34 > 36 > 46 Ar beams at 33 MeV/nucleon 
by a hydrogen target is presented and the dependence on the neutron-proton asymmetry of the 
spectroscopic factors is further investigated. A finite range adiabatic model is used in the analysis 
and compared to previous zero range and local energy approximations. Full three-body Faddeev 
calculations are performed to estimate the error in the reaction theory. In addition, errors from the 
optical potentials are also evaluated. From our new spectroscopic factors extracted from transfer, 
it is possible to corroborate the neutron-proton asymmetry dependence reported from knockout 
measurements. 

PACS numbers: 24.10. Ht; 24.10.Eq; 25.55.Hp 



I. INTRODUCTION 

Since the origin of the shell model [l|, concepts such 
as orbital ordering, level occupancy and magic numbers 
have been extremely helpful to describe nuclear proper- 
ties. In the last couple of decades, with the advances 
of rare isotope facilities, it has been possible to explore 
shell structure away from stability and perhaps expect- 
edly, new orbital ordering, level occupancies and magic 
numbers have been revealed (e.g. the doubly magic 24 
0). In this quest of understanding shell evolution as 
one moves toward the limits of stability, nuclear reac- 
tions have played a central role. Examples of the ongoing 
intense activity are: the systematic knockout program 
at the National Superconducting Cyclotron Laboratory 

(nscl) nm, the sequence of transfer measurements on 

unstable fission fragments at Oak Ridge National Lab- 
oratory (e.g. the experiments by the High Res- 
olution Array (HiRA) collaboration at NSCL (e.g. 0) 
and the transfer program at Argonne National Labora- 
tory on light systems (e.g. [13) EH)- It is of paramount 
importance that the reaction theory used in the analysis 
of these reactions provide an accurate description of the 
process and that theoretical uncertainties be well under- 
stood. 

Traditionally, the linkage between structure and reac- 
tions has often been made through the notion of a spec- 
troscopic factor, a quantity directly associated with the 
occupancy of the single particle levels [l2j]. By compar- 
ing the theoretical predictions for the cross section to 
the data, experimental spectroscopic factors (S exp ) can 
be extracted. These S exp can then be compared to pre- 
dictions from structure models Sth- Although both, S exp 
and Sth, are model dependent, since the early days of 
the field, the methods for extracting S exp were tuned such 
that S^ xp for adding/removing a nucleon to a closed shell 
nucleus agreed with the independent particle picture. 
This is nicely illustrated in the systematic study pre- 
sented in [13]. This procedure has now been challenged 
by the results from (e, e'p) data [3,[II|]. The (e, e'p) data 



are consistent with a reduction of proton spectroscopic 
factors for many closed shell nuclei by roughly 60%. As 
shown in [15j . the apparent discrepancy between trans- 
fer and electron knockout can be understood by using a 
consistent single particle radial behavior for the removed 
particle as well as an improved reaction theory. Along 
those lines, a re-analysis of (d,p) and (p,d) data on Ca 
isotopes, using Hartree-Fock densities to constrain the 
single particle states, produce spectroscopic factors sim- 
ilar to those from (e, e'p) 

With the possibility to study systems with large pro- 
ton/neutron asymmetry, knockout studies revealed a 
strong dependence of the reduction factor R = S exp / S t h 
on the relative energy between valence neutrons and va- 
lence protons (e.g. [17|). This exciting work has been 
systematized through plotting the reduction factor R as 
a function of the difference between separation energies: 
AS = S p — S n (S n — S p ) for proton (neutron) knock- 
out. In this plot, all the stable systems previously stud- 
ied through (e, e'p) correspond to AS ~ 0, but nuclear 
knockout measurements completed the plot for strongly 
negative and positive values. It was found that when re- 
moving loosely bound particles the reduction factor was 
closer to one, while when removing deeply bound parti- 
cles the reduction factor could be as small as 0.25. 

What is learnt from the analysis of the knockout data 
is that, for removing valence nucleons from closed shell 
nuclei, there is a considerable reduction of the spectro- 
scopic strength as compared to the independent particle 
model, and that most large scale shell model calculations 
account for only part of this reduction. This is a clear 
demonstration of the key role of NN correlations within 
the nucleus. The exact nature of these correlations is 
still not fully understood as well as its dependence on 
proton/neutron asymmetry (see (TR [l9j for two recent 
contributions on this topic). 

Motivated by the exciting results from nuclear knock- 
out experiments, transfer reactions 34 ' 36 ' 46 Ar(p,d) were 
performed at 33 MeV/u, with the intention of studying 
the dependence of the reduction factor on the asymme- 



2 



try R(AS). Results published in Q show a clear dis- 
agreement on the trend of the reduction factor: while 
in knockout experiments one observes a strong variation 
of the reduction factor, in transfer data no such depen- 
dence is present. It has been pointed out that different 
beam energies could be at the heart of the disagreement 
since nuclear knockout experiments were all performed at 
around 70 MeV/u while the transfer measurement were 
done at 33 MeV/u. Consequently, the 34,36,46 Ar(p,d) ex- 
periments have been repeated at higher energy and the 
analysis is in progress. Despite this aspect, from the 
theoretical point of view, it is pertinent to inspect the 
reaction theory used in the analyzes and evaluate the 
uncertainties. 

Due to the renewed experimental interest, there have 
been a number of studies on reaction theory for (d,p) and 
(p,d), particularly concerned with uncertainties in the 
methods. Examples include the study of optical potential 
uncertainties [20(, the effect of the reaction mechanism 
2\\ , the uncertainties in the single particle state [22|, [23[ , 
the description of the deuteron continuum [24| and the 
inclusion of full three-body dynamics in the reaction [25[ . 

As in previous systematic studies [26| . the recent 
34,36,46 Ar(p,d) data [|| were analyzed using the adia- 
batic wave model (ADWA) developed by Johnson and 
Soper [27j]. This model provides an improvement over 
distorted wave Born approximation (DWBA) in that it 
takes deuteron breakup into account explicitly, however 
it also makes a zero-range approximation, reducing the 
accuracy of the method. In [9|, a correction to the zero- 
range approximation is introduced, namely the local en- 
ergy approximation (LEA) [28|. Note that, although not 
as simple as the zero-range ADWA, there is a finite range 
formulation of ADWA [29|, [3(| which has recently been 
studied in detail [3l| and compared to its zero range coun- 
terpart. In [3l| finite range effects are shown to be im- 
portant for (d,p) reaction at deuteron energies around 50 
MeV or larger and the LEA prescription is shown to be 
inaccurate for energies well above the Coulomb barrier. 
While the adiabatic model used in Q is very appealing 
due to its simplicity and the possibility of extracting a 
single spectroscopic factor, one should remember that it 
presents an approximate formulation to the full three- 
body problem p + n + A. The ADWA formulation is 
best if only small n-p distances contribute to the trans- 
fer cross section. Faddeev techniques used to solve the 
exact three-body problem are useful to estimate the un- 
certainty involved in the practical ADWA methods. 

In this work we will re-analyze the results from Q 
with the finite-range version of ADWA and using this 
improved reaction theory obtain the reduction factor for 
34 ' 36,46 Ar (Section II). We will also estimate uncertain- 
ties associated with the reaction theory, namely in what 
pertains the interaction and in the three-body reaction 
dynamics (Section III). A final discussion and conclusions 
will be presented in Section IV. 



II. EXTRACTING SPECTROSCOPIC FACTORS 
FROM THE AR(P,D) DATA 

A. Theory 

The most recent analysis on A(p,d)B reactions [l3l.lT6l. 
HI, HH l33ll have relied on the Johnson and Soper reaction 
model [27]. Given the relative importance of deuteron 
breakup in (d,p) and (p,d) reactions, Johnson and Soper 
consider the three-body problem of n + p + B to describe 
the process. The exact T-matrix for the process can be 
written in prior form as (T^ |: 



T = + U pB - U pA \<t> nB x { +2) > 



(1) 



where <p n B is the initial bound state of the nucleus of in- 
terest A, x p + a describes the relative motion of the proton 
and the target in the initial channel distorted by an opti- 
cal potential U p a, and <f(~) is the exact three-body wave 
function. Johnson and Soper expand the exact three- 
body wave function in a complete set of eigenstates of 
the Hamiltonian of the n + p subsystem: 

(r, R) = Mr)Xd(R) + I (r)Xk(R) , (2) 



which includes the ground state of the deuteron <fid{r) 
and all scattering states (f>k(r) of the n+p system. The 
Jacobi coordinates (r , R) are the vector connecting the 
neutron and the proton, and the center of mass of nu- 
cleus B and the center of mass of the deuteron, re- 
spectively. In Ref.[27], the zero-range approximation 
is made V np (r)(j)d{r) = Dq5(t). With this approxima- 
tion, the three-body problem for calculating \]/( + )(r, R) 
reduces to the solution of an optical-model-like equation 
where the distorting potential is the sum of the nucleon- 
target potentials evaluated at half the deuteron energy 
U Z r{R) =U nA (R) + U pA (R): 



(T R + U ZR (R) -E- e d )xd(R) = 



(3) 



Here Tr is the kinetic energy operator associated with 
R. This distorting potential Uzr incorporates deuteron 
breakup effects within the range of V np and thus is not 
meant to describe deuteron elastic scattering. 

In the Johnson and Soper model, the remnant term 
(U P B — U p a) in EqfT]is neglected and the transition am- 
plitude simplifies to: 



T 



D {xl R{ -\RMnB[R)x p + 2{R')) 



(4) 



;R- 



where R' = 7 — ^~ , 

(m B +m„) 

It is well understood that generally the zero-range ap- 
proximation for the deuteron is not accurate. A full fi- 
nite range version of the adiabatic model of [27| was in- 
troduced by Johnson and Tandy [3(3]. Therein the full 
three-body wave function is expanded in Weinberg states 
of the n-\-p system, and the Schrodinger equation for the 
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three-body problem is reduced to a set of coupled chan- 
nel equations. The full coupled channel equations can 
be solved exactly [34( but simplify tremendously when 
only the first term of the Weinberg expansion is required. 
Nguyen et al. (3lj use this approach to perform a system- 
atic study of the effects of finite range, involving a wide 
range of targets and beam energies. Those results show 
that finite range effects can be very significant at interme- 
diate energies. A popular finite-range correction, the so- 
called local energy approximation (LEA), was introduced 
by Buttle and Goldfarb [28] and relies on the truncation 
to first order of a series dependent on incoming momen- 
tum. Comparisons performed in (3lJ show that the local 
energy approximation performs well for sub-Coulomb en- 
ergies but may fail at intermediate energies. 

Here we follow the procedure in [3l[ and evaluate the 
deuteron distorting potential, taking into account the fi- 
nite range of the deuteron: 

UFR{rl) = — . (5 J 

{<t>d{r)\V np \<p cl { i r)) 

Once the adiabatic deuteron wave XdB is calculated using 
Ufr(R), the transfer amplitude can be obtained from: 
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(+)\ 



((/>dXdB\Vnp\(l>nBXpA 



(6) 



where the remnant term has been neglected (U p b — U p a)- 
Remnant term contributions for the cases under study 
here are no larger than 2%. 



B. Results 

We perform finite range calculations for the three re- 
actions under study: 34 ' 36 ' 46 Ar(p,d) 33 ' 3M5 Ar(g.s.) at 33 
MeV/nucleon as in [3l[. For comparison with Q, we 
include results for the zero-range calculations and the lo- 
cal energy approximation. For the description of the Ar 
bound states of interest, we use a Woods Saxon mean 
field with radius r = 1.25 fm and diffuseness a = 0.65 
fm as in [9|, plus a spin-orbit interaction V so — 6 MeV 
with the same geometry as the mean field (note that the 
spin-orbit force in the bound state does not affect the ex- 
tracted spectroscopic factors). Details on the single- hole 
bound states can be found in Table fl] 



TABLE I: Properties of the single hole states: S„ is the sepa- 
ration energy, Rrms is the r.m.s. radius of the hole state, and 
b n ij is it asymptotic normalization constant [22l ]. 



nucleus 


state 


Sn (MeV) Rrms \b„lj\ 


33 Ar 
35 Ar 
45 Ar 


2si/ 2 
W3/2 
I/7/2 


17.07 3.39 31.8 
15.25 3.40 10.3 
8.07 4.10 2.35 
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Angular distributions for: (a) 
33 MeV, (b) 36 Ar(p,d) 35 Ar(g.s.) 



FIG. 1: (Color online) 
34 Ar(p,d) 33 Ar(g.s.) E v = 

E p = 33 MeV and (c) 46 Ar(p,d) 45 Ar(g.s.) E p = 33 MeV. 
Comparison of full finite range (solid) with the zero-range 
approximation (dashed), and the local energy approximation 
(long-dashed). All distributions have been multiplied to scale 
the data by the indicated spectroscopic factor S. 
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TABLE II: Extracted spectroscopic factors using the finite- 
range ADWA model SF(ADWA-RF), and the percentage dif- 
ferences with results obtained with the zero-range approxi- 
mation, the local energy approximation, and the inclusion of 
the remnant term. Also shown are the spectroscopic factors 
calculated in large scale shell model SF(LSSM) 



Model 


34 Ar 36 Ar 46 Ar 


SF(ADWA-FR) 
SF(LSSM) 


0.92 2.21 4.93 
1.31 2.10 5.16 


zero range 
local energy approx. 
remnant 


5% 9% 5% 
14% 9% 3% 
2% 0.3% 0.6% 



understand the uncertainties introduced with the several 
approximations made. The most important approxima- 
tion of course is that we transform the many-body prob- 
lem into a three-body problem, which introduces effective 
interactions and with them uncertainties. However, there 
are also approximations at the three-body level, namely 
the truncation of the expansion in Weinberg states as well 
as not taking into account the coupling between breakup 
and transfer channels to all orders. The exact solution 
to the full three-body can be obtained in the Faddeev 
framework. Below we present some of the essential fea- 
tures of the Faddeev method in momentum space and 
follow with a discussion on the estimation of reaction 
theory uncertainties for the problem. 



For the nucleon optical potentials we use Chapel Hill 
[37| as in Q unless otherwise stated. Finally, for the 
deuteron bound state and the V np interaction appearing 
in the transfer operator of Eq.©, we use the Reid inter- 
action [36] and the corresponding D in the zero range 
model. As stated earlier, the remnant term is neglected 
but was determined to be less that 2% in the cases here 
studied. The finite range deuteron adiabatic potentials 
of Eq.© are generated using the code twofnr [39| and 
the transfer cross sections are computed with the code 

FRESCO [40| . 

Full finite range calculations (solid line) are compared 
to the zero-range approximation (dashed line) and the 
local-energy approximation (long-dashed line), together 
with the measured angular distributions in Fig[TJ All 
calculations have been renormalizcd to the data using 
the same procedure as in Q: the first 3 data points are 
used for 34 Ar, the 5 points around the peak are used 
for 36 Ar and all points but the first are used for 46 Ar. 
The extracted spectroscopic factors S are provided in 
the legend of FigJTJ Overall the angular distributions are 
well described by theory. It is only for the 36 Ar case 
that we find finite range effects to introduce differences 
in the shape of the distribution. For this case these ef- 
fects improve the description of the data. Most impor- 
tantly, the normalization of the cross section changes sig- 
nificantly with finite range effects. As compared to cross 
sections calculated including full finite range effects, zero 
range results systematically overestimate the cross sec- 
tion: 5% for 34 Ar(p,d) 33 Ar, 9% for 36 Ar(p,d) 35 Ar and 
5% for 46 Ar(p,d) 45 Ar. On the other hand, LEA results 
underestimate the cross section by 14% for 34 Ar(p,d) 33 Ar 
and 9% for 36 Ar(p,d) 35 Ar while overestimating the cross 
section by 3% for 46 Ar(p,d) 45 Ar. 



III. ESTIMATING ERRORS FROM REACTION 
THEORY 

The finite range adiabatic model is by no means the 
full solution and an important part of the analysis is to 



A. Faddeev theory 

Recently, nuclear reactions involving protons or 
deuterons have been studied solving the full Faddeev- 
type equations in momentum space |4lj |. These equa- 
tions, usually referred to as AGS equations (for Alt, 
Grassberger and Sandhas) [42| , are coupled integral equa- 
tions for the transition operators 

3 

Upa = (1 - S Pa )G^ + - 6 0a )T a G Q U aa (7) 

where a, f3 refer to the three Faddeev components, which 
also denotes the associated two-body pairs. The on-shell 
matrix elements (ipp\Up a \ip a ) are scattering amplitudes 
and therefore lead directly to the observables. In Eq. ([7]) 
Go = (E + iO — i?o) 1 is the free resolvent, E being the 
available three-particle energy and Hq the free Hamilto- 
nian. The two-particle transition operators T a are ob- 
tained from the Lippmann-Schwinger equation 

T a = v a + v a G T a (8) 

where v a is the potential for the pair a in the odd- 
man-out notation. Thus, AGS equations may be viewed 
as a way for summing up the multiple scattering series 
in terms of T a . Each Faddeev channel state \tpp) for 
13 = 1, 2, 3 is an eigenstate of the channel Hamiltonian 
Hp = Ho + vp with the energy eigenvalue E; thus, \ipp) 
is a product of the bound state wave function for pair 
(3 and a plane wave with fixed on-shell momentum cor- 
responding to the relative motion of particle j3 and pair 
f3 in the initial or final state. The channel states |^o) 
are the eigenstates of Hq with the same eigenvalue E 
and describe the free motion of three particles. Observ- 
ables of elastic scattering are calculated from the ma- 
trix elements with /3 = a, those of breakup are given 
by (3 = while ^ j3 ^ a correspond to transfer reac- 
tions. Although the original AGS equations were derived 
for short-range potentials, the Coulomb interaction can 
be included using the method of screening and renor- 
malization J43|. The numerical details of solving AGS 
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equations with potentials of general form can be found in 
QUI! . Since elastic scattering, transfer, and breakup 
arc treated on equal footing, the Faddeev/AGS frame- 
work yields the most accurate solution to the general 
three-body problem where all channels (elastic, transfer 
and breakup) are fully coupled; it already has provided 
important tests to other methods commonly used in re- 
action theory [25|, |46|, |47[ ■ 



B. Discussion on uncertainties from reaction 
theory 

Although the Faddeev/AGS framework offers an ex- 
act solution for elastic scattering, transfer and breakup, 
given a fixed three-body Hamiltonian, it does have sev- 
eral conceptual limitations regarding applications to nu- 
clear reactions because these are not genuine three-body 
problems. A first one is related to the mapping onto the 
many-body problem. As pointed out in |35| , one of the 
important features of the transition amplitude for the 
transfer process where only V np appears in the transfer 
operator is that the mapping onto the many-body prob- 
lem can be done in a rather straight-forward manner 
in terms of one-nucleon overlap functions. In the Fad- 
deev approach this is no longer possible because elastic, 
breakup and transfer components are all mixed. This im- 
plies that the Faddeev approach to (p,d) reactions should 
only be used for pure single particle states of nucleus A 
and cannot be used to extract a spectroscopic factor from 
a ratio to experimental cross sections (5lj . In this work 
we will use the Faddeev solutions to estimate the theo- 
retical errors in the finite range adiabatic model only. 

Second, the rigorous Faddeev theory can only use the 
same potentials in the initial and final states, while in the 
adiabatic approaches for A(p,d)B reactions [13, HH 
the nuclcon-nuclcus interaction Vnb is different in the 
initial and final channel. Within the adiabatic ap- 
proaches, for the neutron in the initial channel we use 
an effective real binding potential fitted to the spectrum 
of nucleus A whereas in the final state the nucleon opti- 
cal potential contains absorption which reproduces elas- 
tic scattering. In previous Faddeev A(p,d)B calculations 
V n B was real in all partial waves [H, |4l| while V p b was 
complex with parameters corresponding to the A + p 
state. In the present work, Faddeev calculations take 
V n B to be real in the partial wave with the n 4- B bound 
state but in all other partial waves V n B is taken to be 
complex as in the adiabatic approach in the final state. 
The effect of neutron absorption in the partial wave with 
the n+B bound state can be estimated by introducin g en - 
ergy dependent interactions into the AGS equations [48| . 
Even if not fully consistent, we include these calculations 
here to provide an estimate of the magnitude of the ef- 
fect of energy dependence in the neutron- Ar interaction. 
Within adiabatic approaches, the proton interaction in 
the initial channel U p a is calculated at the beam energy, 
whereas in the exit channel U p b it is calculated at half 



the deuteron energy. In Faddeev calculations, the proton 
interaction U p b is fixed at the beam energy. 

There are also some technical challenges in the Fad- 
deev/AGS framework regarding applications to nuclear 
reactions, related to the convergence of multiple scatter- 
ing series, partial wave expansion and Coulomb screen- 
ing. However, these difficulties do not preclude us from 
obtaining well converged results for the reactions consid- 
ered in the present work. 

Finally, another important consideration in determin- 
ing the uncertainties in reaction theory relates to the 
optical potentials. The finite range adiabatic theory is 
built on nucleon optical potentials which are much bet- 
ter known than deuteron optical potentials. Nevertheless 
these are not uniquely determined and thus, in addition 
to the reference Chapel Hill potential 37[, we also per- 
formed calculations with a more recent global potential 
by Koning and Delaroche [Hj]. From the percentage dif- 
ference in cross sections obtained with these two global 
potentials, we estimate the optical potential uncertainty 
in the theoretical cross section. 



C. Estimates of the uncertainties from reaction 
theory 

In order to estimate the error from reaction theory, 
full three-body Faddeev calculations were performed us- 
ing the same interactions, as discussed in Section III Bl 
and IIIIBI The results from these calculations (FADD) 
are compared to the finite range adiabatic calculations 
(ADWA-FR) in Fig. [2] The solid green lines correspond 
to the Faddeev results, while the solid black lines corre- 
spond to the finite range adiabatic. In these plots, no 
normalization to the data is performed. Percentage dif- 
ferences between cross sections obtained within the finite 
range adiabatic model and the Faddeev equations are 
calculated relative to the Faddeev results. These com- 
parisons are performed in the angular region where spec- 
troscopic factors were extracted from the data Q and are 
used to estimate the error of the approximations in the 
reaction mechanism. The relevant angles are: 9 « 9° for 
34 Ar, 6 » 20° for 36 Ar and 6 w 5° for 46 Ar. The percent- 
age differences are presented in Table HTH In addition, we 
also perform Faddeev calculations where instead of using 
the Chapel Hill nucleon potential, we use the more mod- 
ern Koning and de la Roche 



The resulting angular 
distributions (FADD-KD) correspond to the dashed red 
lines in Figs. HJ Percentage differences between these 
two Faddeev results in the same angular range of inter- 
est are presented in Table IIIII For completeness, and to 
ensure that our error estimates in Table IIIII are reliable, 
we also perform Faddeev calculations including the imag- 
inary term of V n B also in the partial wave with the n + B 
bound state, for all E nB > [48]. Results (FADD-E) are 
shown in Fig. [5] as the blue dashed lines and we find that 
differences between these two Faddeev results are of the 
same magnitude as the errors estimated when the V n B is 
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FIG. 2: (Color online) Angular 

E p = 33 MeV, (b) 
46 Ar(p,d) 45 Ar(g.s 



34 Ar(p,d) 33 Ar(g.s 
E p = 33 MeV and (c) 



distributions for: (a) 
36 Ar(p,d) 35 Ar(g.s.) 



E,, 



33 MeV. 



Full three-body Faddeev calculations (FADD) with the finite 
range adiabatic model (ADWA-FR). More detail in the text. 



real only. Since there are ambiguities when introducing 
the energy dependence of V n B, we do not use FADD-E 
to estimate errors. 



TABLE III: Estimates of theoretical errors in the extracted 
spectroscopic factors due to approximations in the reaction 
model as well as experimental errors. 



Errors 


6 th ( 34 Ar) e tfc ( 36 Ar) e th ( 46 Av) 


Optical potential 
Faddeev 
Experiment 


8 % 7% 4% 
6 % 19% 11% 
8% 8% 8% 


Total 


13 % 22 % 14 % 



As the optical potential uncertainty is independent of 
the three-body effects included in the Faddeev descrip- 
tion, and these are both independent from the experi- 
mental statistical error, we estimate the total error by 
summing in quadrature. These are shown in the last line 
of Table IIIII The largest increase in error is seen for the 
36 Ar case, and originates from the approximations in the 
reaction mechanism. 

Finally, as in previous studies Q, we calculate the re- 
duction factor R s as a ratio of the experimentally deter- 
mined spectroscopic factor and that obtained from large 
basis shell model calculations. For the purpose of consis- 
tency, we use the same shell model calculations reported 
in 0. These include the sd-shell model space and the 
USDB effective interaction [49( for 34 > 36 Ar and the sd- 
pf m odel space with the interaction of Nummela et al. 
50] for 46 Ar. The predicted ground state spectroscopic 
factors for 34 Ar, 36 Ar and 46 Ar are 1.31, 2.10 and 5.16, 
respectively, as shown in Table II. The reduction factors 
for 34 Ar, 36 Ar and 46 Ar using the spectroscopic factors 
extracted from transfer cross sections and finite range 
adiabatic theory are plotted as black squares in Fig. [3] 
and errors bars include only the experimental error. The 
total errors, including those from Table IIIII are repre- 
sented by the green bars. The knockout results [|| are 
plotted as red circles and only include statistical errors. 



IV. CONCLUSIONS 

In summary, we have reanalyzed the reactions 
34,36,46 Ar(p,d) 33 ' 35 ' 45 Ar at 33 MeV/nucleon using an im- 
proved reaction theory which includes deuteron breakup 
and finite range effects, based on the adiabatic wave ap- 
proximation (ADWA). In addition we have quantified the 
errors in the reaction theory due to the optical potential 
and the approximate solution of the three-body problem. 
In order to do this we have performed exact three-body 
Faddeev calculations. Since our theory is based on nu- 
cleon optical potentials, we find that the optical potential 
uncertainty is below 10%. However, differences between 
the adiabatic results and the Faddeev results can be as 
large as ss 20%. Further detailed studies to better un- 
derstand the source of these differences are underway. 

With the present level of accuracy of the reaction the- 
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FIG. 3: (Color online) Reduction factors Rs = SF(ADWA — 
FR)/SF(LB — SM) as a function of the difference between 
the neutron and proton separation energies AS. The squares 
and circles correspond to values extracted using transfer or 
knockout respectively. The bars correspond to the total un- 
certainty including both experimental and theoretical errors 
evaluated for the transfer reactions. 



ory, the slope suggested by the knockout data can be 



corroborated by the transfer data. Our error bars do 
not include the uncertainties in the single neutron-hole 
overlap functions. Here we use standard geometry and 
it is not clear whether, specially for 4 6Ar, this will be 
an adequate choice. Measurements of the matter radius 
could provide additional constrains. Ideally, these reac- 
tions could be repeated at lower energy in order to extract 
the asymptotic normalization of the overlap functions, 
which would then enable a much better handle on this 
additional ambiguity [22| . 

Although in this work we concentrate on the reaction 
theory for transfer, it is equally important to estimate 
theoretical errors associated with the description of the 
knockout reactions. Our results call for a better under- 
standing of the reaction mechanism in order to reduce 
the errors. 
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